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Abstract —We consider three samples of O- and B-type stars from the solar 
neighborhood 0.6-4 kpc for which we have taken the distances, line-of-sight velocities, 
and proper motions from published sources. The first sample contains 120 massive 
spectroscopic binaries. O stars with spectroscopic distances from Patriarchi et al. 
constitute the second sample. The third sample consists of 168 OB3 stars whose distances 
have been determined from interstellar calcium lines. The angular velocity of Galactic 
rotation at the solar distance fig, its two derivatives fig and fig, and the peculiar velocity 
components of the Sun (U, V, W) & are shown to be well determined from all three samples 
of stars. They are determined with the smallest errors from the sample of spectroscopic 
binary stars and the sample of stars with the calcium distance scale. The fine structure of 
the velocity held associated with the influence of the Galactic spiral density wave clearly 
manifests itself in the radial velocities of spectroscopic binary stars and in the sample of 
stars with the calcium distance scale. 


INTRODUCTION 

Young massive stars of spectral types O and B are visible from great distances from the 
Sun. Therefore, they are an important tool for studying the structure and kinematics of 
the Galaxy (Moffat et al. 1998; Zabolotskikh et al. 2002; Avedisova 2005; Popova and 
Loktin 2005). Since OB stars do not recede far from their birthplaces in their lifetime, 
they trace well the spiral structure of the Galaxy (Efremov 2011; Coleiro and Chaty 2013; 
Vallee 2014; Hou and Han 2014). 

So far the distances to active star-forming regions, where OB stars are concentrated, 
have been determined by the kinematic method (Anderson et al. 2012), which has a low 
accuracy (Moises 2011). The spectroscopic and photometric distance estimation methods 
can be used only for a relatively small solar neighborhood with a radius of 4-5 kpc. 
However, several factors presently reduce the accuracy of these methods. First, there 
are virtually no direct measurements of the parallaxes for supergiants in the Hipparcos 
catalogue (1997) with an accuracy of at least 10%. Such measurements are needed to 
establish a reliable calibration. Second, the uncertainty in establishing the spectral type 
and luminosity class is great for high-luminosity stars. According to Wegner (2007), the 
first two factors cause, for example, the absolute magnitude dispersion for the sequences of 
la or lab supergiants on the Hertzsprung-Russell diagram to be about 1.5-2.0 magnitudes 
(this dispersion is only about 0.3 magnitude for main-sequence stars). Third, for the 
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application of photometric methods (for example, based on Strorngren photometry), there 
are simply no data for distance determinations. 

The appearance of infrared observations, in particular, the 2MASS catalogue 
(Skrutskie et al. 2006), has contributed to an increase in the accuracy of these 
methods when the interstellar extinction is taken into account. Using the photometric 
characteristics from this catalogue, Patriarchi et al. (2003) determined the spectroscopic 
distances for 184 O stars. One of our goals is to study the Galactic kinematics using data 
on these stars. 

Spectroscopic binaries are of interest in their own right. They usually have a long 
history of observations. The systemic line-of-sight velocities, spectral classification, and 
photometry are well known for many of them. Having analyzed the detached spectroscopic 
binaries with known spectroscopic orbits selected on condition that the masses and radii 
of both components were determined with errors of no more than 3%, Torres et al. (2010) 
showed their spectroscopic distances to agree with the trigonometric ones within error 
limits of no more than 10%. In Eker et al. (2014), the number of such binaries was 
increased by 67%; it is 257. However, the fraction of young massive binaries among them 
is small (about 20 binaries). 

Previously (Bobylev and Bajkova 2013a), we presented a kinematic database 
containing information about 220 massive young stars. The stars are no farther than 
3-4 kpc from the Sun; about 100 of them are spectroscopic binary and multiple systems 
whose components are massive OB stars, while the rest are B stars from the Hipparcos 
catalogue with spectral types no later than B2.5 and errors in the trigonometric parallaxes 
no more than 10%. All stars have estimates of their proper motions, line-of-sight velocities, 
and distances, which allows their space velocities to be analyzed. A significant fraction 
of the stars from this sample (162 stars) belongs to the Gould Belt (with heliocentric 
distances r < 0.6 kpc). Previously (Bobylev and Bajkova 2013a, 2014b), we analyzed the 
kinematics and spatial distribution of these nearby stars. In contrast, the number of more 
distant stars (58 stars at r > 0.6 kpc) should be increased for a reliable determination of 
the parameters describing the Galactic kinematics. 

At present, there are estimates of the distances to OB stars made by the spectroscopic 
method from the broadening of interstellar Ca II, Na I, or K I absorption lines (Megier et 
al. 2005, 2009), with this scale having been reconciled with the trigonometric parallaxes 
(Megier et al. 2005). Quite recently, new measurements of stars performed by this method 
have been published (Galazutdinov et al. 2015). 

The goal of this paper is to expand the kinematic database on massive young star 
systems, to compare the various distance scales, and to determine the parameters of the 
Galactic rotation and spiral density wave based on the data obtained. 

METHOD 

From observations we know three projections of the stellar velocity: the line-of-sight 
velocity V r as well as the two velocities V) = 4.74r/pcos& and Vj, = 4.74r /4 directed along 
the Galactic longitude / and latitude b and expressed in km s' 1 . Here, the coefficient 4.74 
is the ratio of the number of kilometers in an astronomical unit to the number of seconds 
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in a tropical year, and r is the heliocentric distance of the star r in kpc. The proper 
motion components yicosb and fib are expressed in milliarcseconds per year (mas yr -1 ). 
The velocities U, V, and W directed along the rectangular Galactic coordinate axes are 
calculated via the components V r ,Vi, and 14 : 

U — V r cos l cos b — Vi sin l — V, cos l sin b , 

V — V r sin / cos b + Vi cos l — 14 sin l sin 6 , ( 1 ) 

W — V r sin 6 + 14 cos b , 

where U is directed from the Sun to the Galactic center, V is in the direction of Galactic 
rotation, and W is directed toward the north Galactic pole. We can find two velocities, 
Vr directed radially away from the Galactic center and V c i rc orthogonal to it and pointing 
in the direction of Galactic rotation, based on the following relations: 

Vcirc = u sin 9 + (V 0 + V) cos 9, , , 

V R = —U cos 9 + (Go + V) sin 9, 1 j 

where the position angle 9 is calculated as tan 9 = y/(Ro — x ), while x and y are the 
rectangular Galactic coordinates of the star, flo is the angular velocity of Galactic rotation 
at the solar distance Ro , the parameters Oq and il" are the corresponding derivatives of 
the angular velocity, and V 0 = |i? 0 ^o|- 

The velocities Vr and W are virtually independent of the pattern of the Galactic 
rotation curve. However, to analyze any periodicities in the tangential velocities, it is 
necessary to determine a smoothed Galactic rotation curve and to form the residual 
velocities A V C i rc - As experience shows, to construct a smooth Galactic rotation curve in 
a wide range of distances R , it is usually sufficient to know two derivatives of the angular 
velocity, Hq and fig. Note that all three velocities Vr , A V C i rc , and W must be freed from 
the peculiar solar velocity U & ,V Q , W & . 

To determine the parameters of the Galactic rotation curve, we use the equations 
derived from Bottlinger’s formulas in which the angular velocity fl was expanded in a 
series to terms of the second order of smallness in r/R 0 : 

V r = — U & cos 6 cos/ — V Q cosfesin/ , , 

—W Q sin b + R 0 (R — R 0 ) sin l cos bQ' 0 + 0.5Rq(R — Rq) 2 sin l cos 6 flg, 


V = Uq sin l — V Q cos l — rfl 0 cos b . 

+ (R — R 0 )(R 0 cos/ — r cos 6 )flg + 0.5(1? — i? 0 ) 2 (i ? 0 cos l — r cos 6 )Oq, 

Vb — U e cos l sin b + V e sin l sin b , 

—W Q cos b — Ro(R — Rq) sin / sin bQ' 0 — 0.5R o (R — Rq) 2 sin l sin 60 q, 

where R is the distance from the star to the Galactic rotation axis, 

R 2 = r 2 cos 2 b — 2R 0 r cos b cos l + R 2 0 . (6) 


The influence of the spiral density wave in the radial, Vr, and residual tangential, A V C i rc , 
velocities is periodic with an amplitude of ~10 km s -1 . According to the linear theory of 
density waves (Lin and Shu 1964), it is described by the following relations: 


Vr = ~f R COS X, 

A V circ = fe sin x, 


( 7 ) 
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where 


X = m[cot(i)ln(R/R 0 ) -6] + x& ( 8 ) 

is the phase of the spiral density wave (m is the number of spiral arms, i is the pitch angle 
of the spiral pattern, x& is the radial phase of the Sun in the spiral density wave); }r and 
fg are the radial and tangential velocity perturbation amplitudes, which are assumed to 
be positive. 

In the next step, we apply a spectral analysis to study the periodicities in the velocities 
Vr and AV C i rc . The wavelength A (the distance between adjacent spiral arm segments 
measured along the radial direction) is calculated from the relation 


2ttRq 

A 


m cot (i). 
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Let there be a series of measured velocities Vr u (these can be both radial, Vr and residual 
tangential, A Vg, velocities), n — 1,..., N, where N is the number of objects. The objective 
of our spectral analysis is to extract a periodicity from the data series in accordance with 
the adopted model describing a spiral density wave with parameters Jr, A(i) and Xo- 
Having taken into account the logarithmic character of the spiral density wave and 
the position angles of the objects 6 n , our spectral (periodogram) analysis of the series 
of velocity perturbations is reduced to calculating the square of the amplitude (power 
spectrum) of the standard Fourier transform (Bajkova and Bobylev 2012): 



N 

Y V n( R 'n) eX P 


n=l 
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where V\ k is the kth harmonic of the Fourier transform with wavelength A*, = D/k,D is 
the period of the series being analyzed, 


R’ n = R 0 \n(Rn/R 0 ), 

Vn(R’n) = V n(R’n ) X exp (jm9 n ). 


The algorithm of searching for periodicities modified to properly determine not only the 
wavelength but also the amplitude of the perturbations is described in detail in Bajkova 
and Bobylev (2012). 


DATA 

The Sample of Spectroscopic Binary OB Stars 

To select young massive distant spectroscopic binaries, we used the SB9 database 
(Pourbaix et al. 2004), which contains references to the spectroscopic orbit determinations 
up until 2014. We considered binaries with spectral types of the primary component no 
later than B2.5 and various supergiants with luminosity classes la and lab. Thus, we are 
interested in stars with masses of more than 10 solar masses. 

For a preliminary distance estimate, we consulted various sources, for example, 
the compiled database by Gudennavar et al. (2012) and the FUSE (Far Ultraviolet 


4 



Ta6jiHU,a 1: Bibliographic information about the stars 


Star 

Vy, /x, dist 

Star 

Vj, /i, dist 

MWC 1 

Sb9, L07, G12 

Herschel 36 

ArlO, UC4, ArlO 

V745 Cas 

Chl4a, L07, Chl4a 

V411 Ser 

Sb9, L07, G12 

BM Cas 

Sb9, L07, Pop77 

mu. Sgr 

Sb9, L07, BB15 

HD 12323 

Sb9, UC4, Gall5 

16 Sgr 

Mal4, L07, G12 

V622 Per 

Mal07, UC4, Cap02 

MY Ser 

Ibl3, L07, Ibl3 

V615 Cas 

MJ14, MJ14, MJ14 

QR Ser 

S09, L07, S09 

V482 Cas 

Sb9, Tyc2, Mc03 

RY Set 

Sb9, L07, Sm02 

MY Cam 

Lrl4, Tyc2, Lrl4 

MCW 715 

Sb9, L07, Hut71 

KS Per 

Sb9, L07, BB15 

V337 Aql 

Tul4, Tul4, Tu 14 

HD 31617 

Sb9, L07, Guil2 

MWC 314 

Lol3, Tyc2, CarlO 

HD 31894 

Sb9, L07, Guil2 

VI765 Cyg 

Sb9, L07, G12 

eps Aur 

Sb9, L07, Guil2 

V380 Cyg 

Sb9, L07, Tkl4 

ADS 4072AB 

Sb9, L07, Pat03 

HD 124314 

SIMB, L07, G12 

3 Pup 

Sb9, L07, RR94 

V2107 Cyg 

B14, L07, B14 

HD 72754 

Sb9, L07, L07 

ALS 10885 

G06, L07, Pat03 

ALS 1803 

Sb9, L07, Pat03 

V470 Cyg 

Sb9, L07, Coml2 

ALS 1834 

Sb9, Tyc2, Pat03 

WR 140 

Sb9, Dz09, Do05 

V560 Car 

Sb9, Tyc2, R01 

V404 Cyg 

MJ14, MJ14, MJ14 

V346 Cen 

Sb9, UC4, KalO 

Schulte 73 

Ki09, UC4, Ki09 

KX Vel 

Mal4, L07, Mal4 

V382 Cyg 

Sb9, L07, YY12 

HD 97166 

Mal4, UC4, Mal4 

V2186 Cyg 

Kil2, UC4, Kil2 

HD 115455 

Mal4, UC4, Mal4 

MT91 771 

Kil2, UC4, Ki07 

HD 123590 

Mal4, L07, Av84 

VV Cep 

Sb9, L07, L07 

HD 150136 

S13, UC4, S13 

V2174 Cyg 

Sb9, L07, Bol78 

HD 152234 

Sb9, L07, G12 

4U 2206+54 

StZ14, Tyc2, SIMB 

HD 152246 

Nal4, L07, Nal4 

V446 Cep 

Chl4b, L07, Chl4b 

Braes 144 

SIMB, Tyc2, G12 

MWC 656 

Casl4, L07, Casl4 

V1297 Sco 

SIMB, Tyc2, G12 

GT Cep 

Chl5, L07, Chl5 

Pismis 24-1 

MW07, Tyc2, MW07 

V373 Cas 

Sb9, L07, H187 

HD 163892 

Mal4, L07, Call 

ALS 12775 

Sb9, L07, BB15 

HR 6716 

Sb9, L07, St99 

VZ Cen 

Will3, L07, BB15 


1. The sources of line-of-sight velocities: Sb9 (Pourbaix et al. 2004); Chl4a, Chl4b, Chl5 (Qhakirli et al. 
2014a, 2014b, 2015); ArlO (Arias et al. 2010); Mal4 (Mayer et al. 2014); G06 (Gontcharov 2006); Mal07 
(Malchenko et al. 2007); Ibl3 (Ibanoglu et al. 2013); MJ14 (Miller-Jones 2014); S09, S13 (Sana et al. 
2009, 2013); Lrl4 (Lorenzo et al. 2014); Tul4 (Tiiysiiz et al. 2014); StZ14 (Stoyanov et al. 2014); Ki07, 
Ki09, Kil2 (Kiminki et al. 2007, 2009, 2012); Loi3 (Lobel et al. 2013); B14 (Baki§ et al. 2014); Nal4 
(Nasseri et al. 2014); Casl4 (Casares et al. 2014); MW07 (Maiz-Apellaniz et al. 2007); Will3 (Williams 
et al. 2013); SIMB (SIMBAD database). 

2. The sources of proper motions: L07 (van Leeuwen 2007); UC4 (UCAC4, Zacharias et al. 2012); Tyc2 
(Tycho-2, Hog et al. 2000); Dz09 (Dzib and Rodriguez 2009). 

3. The sources of distances: G12 (Gudennavar et al. 2012); Call (Megier et al. 2009); Pat03 (Patriarchi 
et al. 2003); Pop77 (Popper 1977); Gall5 (Galazutdinov et al. 2015); Cap02 (Capilla and Fabregat 2002); 
Mc03 (McSwain 2003); Sm02 (Smith et al. 2002); Av84 (Avedisova and Kondratenko 1984); Bol78 (Bolton 
and Rogers 1978); Hut71 (Hutchings and Redman 1971); Guil2 (Guinan 2012); CarlO (Carmona et al. 
2010); R01 (Rauw et al. 2001); RR94 (Rovero and Ringuelet 1994); Coml2 (Comeron and Pasquali 
2012); Do05 (Dougherty et al. 2005); H187 (Hill and Fisher 1987); St99 (Stickland and Lloyd 1999); Tkl4 
(Tkachenko et al. 2014); KalO (Kaltcheva and Scorcio 2010); BB15 (our estimates); YY12 (Yagarsoy and 
Yakut 2012). 
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Y. kpc 


Phc. 1: Distribution of 120 spectroscopic binary stars in projection onto the Galactic XY 
plane; the Sun is at the intersection of the dashed lines; the fragments of the four-armed 
spiral pattern are indicated; the position of the local spur, the Orion arm, is marked. 

Spectroscopic Explorer) survey (Bowen et al. 2008). Finally, the publications of various 
authors that determine orbits are an important source of information. 

We did not consider stars nearer than 0.6 kpc, where the Gould Belt dominates, 
and stars farther than 4 kpc, where the proper motion errors are great. We added 62 
more binaries that possess no properties of “runaway” stars (their residual velocities do 
not exceed 40-50 km s” 1 ) to the 58 already described in Bobylev and Bajkova (2013a). 
Bibliographic information about these binaries is given in the table. There are cases where 
there is no orbit in the SB9 catalogue, while it is specified in the SIMBAD database that 
the star is a spectroscopic binary, for example, HD 124314. 

For a number of stars, we calculated the distances by ourselves. The abbreviation BB15 
is used for them in the table. Here, we calculated the distances using the photometric 
magnitudes from the SIMBAD database based on the well-known relation 

logr = 1 + 0.2((V — M v ) — A v ), (12) 

where A v = 3.2 E(B—V), E(B—V) = (B—V) — (B—V) 0 . We took the absolute magnitudes 
of the stars M v from Wegner (2007), where they were determined from Hipparcos stars for 
all spectral types and luminosity classes, and the unreddened colors of the stars (B — V )o 
from Straizys (1977). 

The final sample of spectroscopic binaries contains 120 stars. Their distribution in 
the Galactic XY plane is presented in Fig. 1. The figure shows a fragment of the global 
four-armed (one of the four arms, the Outer one, is not visible here) spiral pattern with 
a pitch angle of 13°. The pattern was constructed in our previous paper (Bobylev and 
Bajkova 2014a), where we analyzed the distribution of Galactic masers with measured 
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H, kpc 


Piic. 2 : (a) Galactic rotation curve constructed from the sample of 120 spectroscopic 

binary stars with parameters (13) with an indication of the boundaries of the la confidence 
intervals; residual rotation velocities A V c i rc (b), radial velocities Vr (c), and vertical 
velocities W (d) of stars versus R ; the vertical dashed line marks the Sun’s position. 
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Pnc. 3: Power spectra of the radial, Vr, residual tangential, AV C i rc , and vertical, W, 
velocities for the sample of spectroscopic binary stars. 


trigonometric parallaxes. The line segment between the Perseus and Carina-Sagittarius 
Arms indicates the Local Arm model (Bobylev and Bajkova 2014b). As can be seen 
from the figure, the distribution of stars agrees well with the plotted spiral pattern. 
Unfortunately, there are few stars in the third quadrant in the Perseus Arm. In this 
paper, we take the Galactocentric distance of the Sun to be Ro = 8.0 ± 0.4 kpc. 

O Stars from the List by Patriarchi et al. (2003) 

Patriarchi et al. (2003) determined the spectroscopic distances for 184 O-type stars; 
photometric data from the 2MASS catalogue were used to take into account the 
interstellar extinction. This sample consists mostly of single stars, although there are 
also spectroscopic binaries in it. 

The proper motions from the Hipparcos catalogue revised by van Leeuwen (2007) or 
from Tycho-2 (Hog et al. 2000) and line-of-sight velocities (according to the SIMBAD 
database) are known for 156 stars from this list. However, the information is complete 
only for 101 of them, i.e., their distances, proper motions, and line-of-sight velocities are 
known simultaneously. The distribution of 156 O stars in projection onto the Galactic 
XY plane is given in Fig. 4. Just as in Fig. 1, the spiral arms are plotted. The agreement 
of the stars from this sample with the plotted spiral pattern is noticeably poorer than 
that for the stars in Fig. 1. Here, we see a larger dispersion of the coordinates in the 
Carina-Sagittarius Arm toward the Galactic center (Z = 0°), in the Perseus Arm toward 
the Galactic anticenter (/ = 180°), and, finally, the Orion Arm segment is noticeably 
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Puc. 4: Distribution of 156 O stars from the list by Patriarchi et al. (2003) in projection 
onto the Galactic XY plane; the Sun is at the intersection of the dashed lines; the 
fragments of the four-armed spiral pattern are indicated; the position of the local spur, 
the Orion Arm, is marked. 



stretched in a direction near l = 75°. 

The Distances Determined from Interstellar Ca II Lines 

The distances to OB stars determined by the spectroscopic method from the broadening 
of interstellar Ca II, Na I, or K I absorption lines are of indubitable interest. The method 
of determining such distances is well known. However, the “calcium” distance scale has 
been tied to the Hipparcos trigonometric parallaxes only recently (Megier et al. 2005). 
The assumption about a uniform distribution of ionized atoms in the Galactic plane 
underlies the method. In some places (for example, toward the cluster Trumpler 16), 
the nonuniformities in the distribution of matter can be significant, which can increase 
the dispersion of the distance estimates by this method. On average, according to the 
estimates by Megier et al. (2009), the accuracy of determining the individual distances to 
OB stars is rH5 %. 

The first catalogue (Megier et al. 2009) contains 290 OB stars, while Galazutdinov 
et al. (2015) determined the distances to 61 more OB stars by this method, with the 
line-if-sight velocities having been measured for all of them. Since the two samples have 
overlaps, the total number of stars in this distance scale is about 340. 

Previously (Bobylev and Bajkova 2011, 2013b), we studied a sample of 258 Hipparcos 
O B3 stars with distances from Megier et al. (2009). About 20% of the sample are 
either known runaway stars or candidates for runaway stars, because they have large 
(> 40 km s -1 ) residual space velocities. As a result, the range of distances 0.8-3.5 kpc 
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Puc. 5: Galactic rotation curve constructed from the sample of 101 O stars from the list 
by Patriarchi et al. (2003) with parameters (14) with an indication of the boundaries of 
the la confidence intervals (a); residual rotation velocities A V c i rc (b), radial velocities Vr 
(c), and vertical velocities W (d) of stars versus R] the vertical dashed line marks the 
Sun’s position. 
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was represented by 102 stars. A slight, no more than 20%, reduction of the distance scale 
was shown to be necessary by analyzing several kinematic parameters obtained from this 
sample (Bobylev and Bajkova 2011). In this paper, we use the entire sample, and not only 
the Hipparcos stars. 

Figure 6 gives the distribution of 168 OB stars with distances in the calcium distance 
scale in the Galactic XY plane for two distance scale factors. In the first case, this factor 
is equal to one; in the second case, the distances to the stars were multiplied by 0.8. Note 
that the new measurements by Galazutdinov et al. (2015) were performed exclusively in 
two narrow sectors: in directions l ~ 135° and l ~ 190° with heliocentric distances of 
about 2.3 kpc. Two structures elongated along the line of sight are clearly seen in Fig. 6 
in these directions. 


RESULTS 

Spectroscopic Binary 

Using 120 spectroscopic binaries, we found the following kinematic parameters from a 
least-squares solution of the system of conditional equations (3)-(5): 

(£/©, V Q , W q ) = (2.8, 9.2, 9.3) ± (1.0, 1 . 2 ,1.0) km s -1 , 

Uo = 29.3 ± 0.8 km s -1 kpc -1 , , . 

Hq = —4.28 ± 0.15 km s -1 kpc -2 , ' ' 

Q" = 0.957 ± 0.128 km s -1 kpc -3 . 

In this solution, the error per unit weight is cr 0 = 10.5 km s -1 . For the adopted i ? 0 = 
8.0 ± 0.4 kpc, the linear Galactic rotation velocity Vo = |i?o^o| is 234 ± 14 km s -1 , while 
the Oort constants are A = —17.1 ± 0.6 km s -1 kpc -1 and B = 12.2 ± 1.0 km s -1 kpc -1 . 
The Galactic rotation curve constructed with parameters (13) and the residual tangential, 
A Vcirc, radial, Vr, and vertical, W, velocities of stars as a function of the distance R are 
given in Fig. 2. 

Why are the vertical velocities W given in Fig. 2? The point is that periodic vertical 
velocity oscillations have been detected quite recently in Galactic masers with measured 
trigonometric parallaxes (Bobylev and Bajkova 2015). Such oscillations with an amplitude 
of about 4-6 km s -1 can be associated with the influence of the Galactic spiral density 
wave. Therefore, confirming this phenomenon in the velocities of other samples of stars 
with different distance scales is of great interest. 

Based on a spectral analysis, we determined the parameters of the spiral density wave 
using the sample of spectroscopic binary stars. The results are reflected in Fig. 3, where 
the power spectra of the radial, Vr, residual tangential, AV%. c , and vertical, W, velocities 
are given. Note that there is a high significance of the signal, p = 0.999, only in the 
power spectrum of the radial velocities. For the tangential velocities, the significance of 
the signal at a wavelength A ~ 2.6 kpc is p = 0.898. As can be seen on the middle panel 
in Fig. 3, there are spurious signals in the range of short wavelengths whose amplitudes 
exceed the amplitude of the signal of interest to us. There are also spurious signals in the 
power spectrum of the vertical velocities W, while the signal at a wavelength of about 1.5 
kpc agrees poorly with the two preceding graphs. 
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Puc. 6 : Distribution of 168 OB stars with distances in the “calcium” distance scale in 
projection onto the Galactic XY plane (a); the distances to these stars were multiplied 
by factor of 0.8 (b); the fragments of the four-armed spiral pattern and the position of 
the Orion Arm are indicated. 


As a result, for the model of a four-armed spiral pattern (m = 4), we found the 
amplitudes of the velocity perturbations Jr = 9.5 ± 1.5 km s -1 and fg = 3.2 ± 1.4 km s -1 
in the radial and tangential directions, respectively, the wavelength Xr = 2.8 ± 0.5 kpc 
(then, iji = —13° ± 4°) and A g = 2.6 ± 0.4 kpc (then, ig = —12° ± 3°) at the Sun’s phase 
in the spiral density wave (%©)/? = —95° ± 15°, and ( Xo)e = —93° ± 12°, respectively. The 
corresponding waves are given on panels (b) and (c) in Fig. 2. Their logarithmic character 
is clearly seen. 

Stars from Patriarchi et al. (2003) 

Using 156 O stars from the list by Patriarchi et al. (2003), we found the following 
parameters: 

([/©, V Q , W q ) = (5.8,16.9, 8 . 6 ) ± ( 2 . 1 , 2.6,1.8) km s -1 , 

O 0 = 29.6 ± 1.2 km s -1 kpc -1 , , . 

Oq = —4.27 ± 0.27 km s -1 kpc -2 , 

Dp = 0.910 ± 0.143 km s -1 kpc -3 . 

In this solution, the error per unit weight is <jq = 21.9 km s -1 , and the linear Galactic 
rotation velocity is Vo = 237 ± 15 km s -1 for R 0 = 8.0 ± 0.4 kpc. 

There are quite a few stars without line-of-sight velocities in this list. Equations (3)- 
(5) can be solved with limited information, which was done when obtaining solution (14). 
However, complete information is needed to calculate the velocities A V c i rc , Vr, and W, 
i.e., only 101 stars can be used. 
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Pnc. 8: Power spectra of the radial, Vr, residual tangential, A V c i rc , and vertical, W, 
velocities for the sample of 0B3 stars with the calcium distance scale. 


The Galactic rotation curve constructed with parameters (14) and the residual 
tangential, A Vd rc , radial, Vr, and vertical, W, velocities of stars as a function of the 
distances R are given in Fig. 5. It can be seen from this figure that a smooth Galactic 
rotation curve is determined quite well. Since the velocities A V c i rc , Vr, and W are fairly 
irregular in pattern, we did not perform a spectral analysis of these velocities. 

The error per unit weight cr 0 characterizes the dispersion of the residuals in the least- 
squares solution of the system of equations (3)-(5). Its value is usually close to the value 
of the residual velocity dispersions for the sample stars averaged over all directions (the 
“cosmic” velocity dispersion). The value of <r 0 = 21.9 km s -1 found in solution (14) is twice 
Co ~ 10 km s -1 expected for OB stars. This may imply that the errors in the distances 
and, consequently, the tangential velocities of these stars are too great. Panel (d) in Fig. 
5, on which large vertical velocities W are seen for a considerable number of stars, attracts 
our attention. 

We produced a new sample of stars with space velocities for two constraints: \W\ <40 
km s -1 and \z\ < 100 pc. Applying these constraints allows the number of possible runaway 
stars to be reduced significantly. 61 stars satisfy these criteria. Using them, we found the 
following parameters: 

(£/©, V Q , W q ) = (10.1,14.6, 5.7) ± (2.1, 2.3,1.9) km s -1 , 

h2 0 = 31.7 ± 1.5 km s -1 kpc -1 , , . 

Uq = —4.53 ± 0.29 km s -1 kpc -2 , 

Oq = 1.124 ± 0.247 km s -1 kpc -3 , 
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where the error per unit weight is <To = 13.7 km s 3 , which is considerably smaller than 
that in solution (14). 

The “Calcium” Distance Scale 

Based on 168 stars, we obtained the following solution for the adopted Ro = 8.0 ±0.4 kpc: 

(U & ,V & ,W Q ) = 

(10.2, 9.3, 9.5) ± (1.1,1.6,1.1) km s -1 , 

= 28.7 ± 0.9 km s -1 kpc -1 , (16) 

fig = —4.21 ± 0.16 km s -1 kpc -2 , 
fig = 0.652 ± 0.149 km s -1 kpc -3 , 

(To = 13.6 km s -1 and Vo = 230 ± 15 km s -1 . Based on the same stars, we also obtained 
the solution for the distance scale factor 0.8: 

(Uq, Vq, W q ) = 

(10.4, 9.7, 7.0) ± (1.0,1.4, 0.9) km s -1 , 

O 0 = 31.3 ± 0.9 km s -1 kpc -1 , (17) 

fig = —4.78 ± 0.16 km s -1 kpc -2 , 
fig = 0.864 ± 0.176 km s -1 kpc -3 , 

(To = 11.6 km s -1 n Vo — 250 ± 15 km s -1 . Figure 7 gives the Galactic rotation curve 
and the velocities A V C i rc , Vr, and W as a function of R constructed from OB stars with 
distances in the original calcium scale and the scale reduced by 20%. 

Based on a spectral analysis, we determined the parameters of the spiral density wave 
from the sample of spectroscopic binary stars with the reduced distance scale. The results 
are reflected in Fig. 8, where the power spectra of the radial, Vr, residual tangential, 
A Vdrc, and vertical, W, velocities are given. Note that there is a high significance of the 
signal, p = 0.999, only in the power spectrum of the radial velocities. It can be clearly seen 
on panels (b) and (c) that there is no statistically significant signal with a wavelength of 
about 2.5 kpc typical of the spiral density wave for the tangential and vertical velocities. 

As a result, from our analysis of the radial velocities we found the velocity perturbation 
amplitude Jr = 11.8 ±1.3 km s -1 and the wavelength A = 2.1 ±0.4 kpc (i = —9.5° ±1.7°) 
at the Sun’s phase in the spiral density wave x© = —86° ±7°. This wave is given on panel 
(g) in Fig. 7. 

DISCUSSION 

First of all, it should be noted that the peculiar solar velocity components, the angular 
velocity of Galactic rotation, and its two derivatives, are determined quite well from all 
three samples of stars. They are determined with the smallest errors from the sample of 
spectroscopic binary stars and the sample of stars with the calcium distance scale. 

For two samples, we also determined such kinematic parameters previously, but only 
using a smaller number of stars. The results of solution (13) should be compared with the 
results of our analysis for 58 distant spectroscopic binary stars outside the circle with a 
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radius of 0.6 kpc (Bobylev and Bajkova 2013a): ( U®,V Q ) = (2.2, 8.4) ± (0.9,1.1) km s -1 
and fl 0 = 31.9 ±1.1 km s -1 kpc -1 , V2g = —4.30 ±0.16 km s -1 kpc -2 , Qg = 1.05 ±0.35 km 
s -1 kpc -3 , where the error per unit weight is cr 0 = 9.9 km s -1 and V 0 = 255 ± 16 km s -1 . 

Based on a sample of 102 OB3 stars with the calcium distance scale (and the scale 
factor 0.8), previously (Bobylev and Bajkova 2011) we found the following parameters: 
(Uq,V & ,W & ) = (8.9,10.3,6.8) ± (0.6,1.0, 0.4) km s -1 , V2 0 = 31.5 ± 0.9 km s -1 kpc -1 , 

= —4.49 ± 0.12 km s -1 kpc -2 , fig = 1.05 ± 0.38 km s -1 kpc -3 , the error per unit 
weight was Oo = 9.6 km s -1 , and the circular rotation velocity was Vo = 252 ± 14 km s -1 
(i?o — 8 kpc). We can see good agreement with solution (17). The only difference is that 
the error in the second derivative of the angular velocity V2 0 decreased considerably with 
increasing number of stars. 

The distances for the O stars from the list by Patriarchi et al. (2003) were determined 
not quite reliably. Note that these authors took the absolute magnitudes of the stars M v 
from the rather old paper by Panagia (1973) and their spectral types from the catalogue 
by Garmany et al. (1982). The correct spectral classification is of great importance in 
determining the spectroscopic distances of stars. A telling example is the star HD 160641, 
for which the spectral type 09.51 is specified in the catalogue by Patriarchi et al. (2003), 
and the distance r — 9 kpc was calculated. According to the estimates by a number of 
authors, this is a helium star, a sdOc9.5II-III:He40 dwarf (Drilling et al. 2013) with a 
mass of ~1 Mq. In this case, the distance to the star is r = 3.3 ± 0.8 kpc (Lynas-Gray et 
al. 1987). 

Why does the sample of spectroscopic binaries look better? After all, the distances 
to the stars were determined either by the photometric method or by the spectroscopic 
one. In our view, the fact that each star is “piece goods” when the spectroscopic orbits 
are determined plays a role here. In this case, the authors have a detailed idea of the 
binary characteristics, in particular, of the spectra of the binary components; therefore, 
the inaccuracies in the spectral classification are minimal. 

At present, there is a sample of approximately 100 maser sources whose trigonometric 
parallaxes have been measured by the VLBI method with a very high accuracy, with a 
mean error of ±20 /zas and, for some of them, with a record error of ±5 /jl as. From an 
analysis of these masers, Reid et al. (2014) found the velocity of the Sun Vo = 240 ± 8 km 
s -1 (i?o = 8.34 ± 0.16 kpc). Previously, based on a smaller number of masers, Honma et 
al. (2012) obtained an estimate of Vo = 238 ± 14 km s -1 (Rq = 8.05 ± 0.45 kpc). It can be 
seen that the values of this velocity found in this paper are in good agreement with the 
best present day estimates. 

Based on a sample of masers, previously (Bobylev and Bajkova 2015) we found the 
following parameters of the spiral density wave for rn = 4 : fo — 6.0 ± 2.6 km s -1 , 
f R = 7.2 ± 2.2 km s - \ \ e = 3.2 ± 0.5 kpc, \ R = 3.0 ± 0.6 kpc, ( X q)o = -79° ± 14° and 
( x&)r — —199° ± 16°. The values found in this paper from the sample of spectroscopic 
binary stars are fairly close to them. The values found from the radial velocities of stars 
with the calcium distance scale are characterized by a larger amplitude and a smaller 
( Xq)r ■ However, this just confirms the results of our previous analysis (Bobylev and 
Bajkova 2011). 

Given the rather significant relative errors in the distances for many of the stars being 
used, the question about the systematic bias in the mean distances of the remaining 
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Pnc. 9: Distribution of 43 stars with the calcium distance scale on the line of sight in the 
direction l = 135°; a Gaussian with a mean of 2.56 pc and a dispersion of 0.54 pc was 
fitted. 


stars arises. Obviously, there are relatively more objects with overestimated distances 
among the rejected stars farther than 4 kpc and, on the contrary, with underestimated 
distances among the stars nearer than 600 pc. Such a bias in the distances is related to 
the Lutz-Kelker effect (Lutz and Kelker 1973; Stepanishchev and Bobylev 2013). For the 
classical case where a uniform distribution of stars is adopted, we can use a formula of 
the form (Lutz and Kelker 1973) 


f(z) = -j ex P 


(z-iy- 




(18) 


where z = vr/7r 0 is the ratio of the true parallax n to the observed parallax 7r 0 , a p = <J 0 /n 0 
is the relative error of the observed parallax, and F(z) is the probability density. 

We took the mean error in the distances for the spectroscopic binaries and the stars 
with the calcium distance scale to be 15%. The mean correction for the Lutz-Kelker effect 
for these stars is then 100 pc (the measured distances should be increased). For the stars 
from the list by Patriarchi et al. (2003), we took the error in the distance for each star to 
be 30%. As a result of our simulations, we found that the total bias for these stars due to 
the Lutz-Kelker effect could be 240 pc. 

The chains of stars elongated along the line of sight, which can be seen in Fig. 4 and 
particularly clearly in Fig. 6, deserve special attention. The most plausible explanation 
is that there are large relative errors in the distances to the objects that actually belong 
to a single compact grouping (a cluster, an association, or a stellar complex). Based on 
such objects, we obtained an independent estimate of the relative error in the distance 
using stars with the calcium distance scale as an example. The result is reflected in Fig. 9, 
where the distribution of 43 OB stars with the calcium distance scale lying on the line of 
sight in the direction l = 135° is given. Based on the parameters of the fitted Gaussian (a 
mean of 2.56 pc and a dispersion of 0.54 pc), we find the relative error in the distance to 
be 21%. The mean value calculated from the individual distance errors for these stars is 
0.42; the relative error in the distance is then 16%. 
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CONCLUSIONS 


We considered 0- and B-type stars located in a wide solar neighborhood, with distances 
from 0.6 to 4 kpc. We produced three samples of stars for which the distances, line-of- 
sight velocities, and proper motions collected from published data were available. The 
first sample includes 120 massive (with masses > 1OM 0 ) spectroscopic binaries. O stars 
for which Patriarchi et ah (2003) estimated the spectroscopic distances using infrared 
photometric data from the 2MASS catalogue constitute the second sample. Single stars 
constitute the overwhelming majority in this sample. The third sample consists of 168 
OB3 stars whose distances were determined from interstellar calcium lines. As in the 
previous case, mostly single stars enter into this sample. 

Such kinematic parameters as the angular velocity of Galactic rotation at the solar 
distance Qg, its two derivatives Qq, and fig and the peculiar solar velocity components 
(U, V, W) & were shown to be well determined from all three samples of stars. They are 
determined with the smallest errors from the sample of spectroscopic binaries and the 
sample of stars with the calcium distance scale. 

The fine structure of the velocity held associated with the influence of the Galactic 
spiral density wave clearly manifests itself in the radial velocities of the spectroscopic 
binaries and in the sample of stars with the calcium distance scale. No small-amplitude 
periodic oscillations in the vertical velocities with a wavelength of about 2.5 kpc manifest 
themselves in any of the stellar samples considered in this paper. The available accuracies 
of the distances are probably insufficient for a reliable detection of such oscillations. 

The linear rotation velocity of the Galaxy Vo was found from the sample of 
spectroscopic binary stars to be 234±14 km s' 1 (for the adopted i? 0 = 8.0±0.4 kpc). Based 
on a spectral analysis of the radial and tangential velocities for these stars, we determined 
the parameters of the spiral density wave. For the model of a four-armed spiral pattern 
(m = 4), we found the following: the perturbation velocity amplitudes fn = 9.5 ± 1.5 km 
s” 1 and fg = 3.2 ± 1.4 km s _1 in the radial and tangential directions, respectively; the 
wavelength A r = 2.8 ± 0.5 kpc (then, in = —13° ± 4°) at the Sun’s phase in the spiral 
density wave (x 0 )r = —95° ± 15° and A g = 2.6 ± 0.4 kpc (then, ig = —12° ± 3°) at the 
Sun’s phase (x&)e = —93° ± 12°. 

The velocity Vo = 250 ± 15 km s ” 1 (for R 0 = 8.0 ± 0.4 kpc and a distance scale 
factor of 0.8) was found from the sample of OB3 stars with the calcium distance scale. 
The parameters of the spiral density wave were determined only by analyzing the radial 
velocities. The perturbation velocity amplitude is /r = 11.8 ± 1.3 km s ' 1 and the 
wavelength is Ar = 2.1 ± 0.3 kpc (m = 4, in = —9.5° ± 1.7°) at the Sun’s phase in 
the spiral density wave (x 0 )r = — 86 ° ± 7°. 
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